Calculation of K ^ tt matrix elements in quenched domain-wall QCD * 
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We explore the possibility for an evaluation of non-leptonic AS = 1 K decay amplitudes through the calculation 
of — > TT matrix elements using domain-wall QCD. The relation between the physical K — > tttt matrix elements 
and K ^ TV matrix elements deduced from chiral perturbation theory is recapitulated. Quenched numerical 
simulations are performed on an 16'^ x 32 x 16 lattice at a lattice spacing ~ 2GeV for both the standard 
plaquette gauge action and a renormalization-group improved gauge action, and reasonable signals for K ~* tt 
matrix elements are obtained. Preliminary results are reported on the K nn matrix elements, and results from 
two actions are compared. 



1. INTRODUCTION 

Despite extensive efforts over the years, lattice 
QCD calculation of the matrix elements relevant 
for understanding the K nn decays, includ- 
ing the AI = 1/2 rule and the direct CP viola- 
tion parameter s'/e^ have achieved only limited 
success to the present [0. One reason behind the 
slow progress is lack of full chiral symmetry in the 
Kogut-Susskind|^ and Wilson|^ formulation for 
lattice fermions employed in the past attempts. 
The domain wall fermion formalism |^,^ offers a 
possibility of resolving this problem. In this arti- 
cle we report on our study of the K tttt ampli- 
tudes in quenched QCD using this formalism^. 

The framework of our study is the reduc- 
tion formulae 01 derived from chiral perturba- 
tion theory (xPT), which relate K ^ t: ampli- 
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tudes to the physical K ^ tttt amplitudes. We 
briefly review these formulae to expose some lim- 
itations |,§. 

Chiral property of domain wall QCD has been 
examined in detail recently jl^,|l^. We have 
shown in particular [p^ that the residual chiral 
symmetry breaking due to finite fifth-dimensional 
size is significantly reduced for a renormalization- 
group improved gluon action as compared to the 
standard plaquette action. We therefore calculate 
the K ^ TT matrix elements for the two gluon ac- 
tions in parallel. 

2. xPT REDUCTION FORMULAE 

The effective Hamiltonian Hw for AS = 1 K 
meson decays is written as 

10 

Hlf^ = ^VudV:^ J2 w{^^hQ^{^^), (i) 

^ i—l 

where W{ii)i are Wilson coefficients. The basis 
four-quark operators Qi are decomposed as Q^^^ + 

(2) 

Ql corresponding to iso-spin in the final state 



2 



/ = and 2. 

Under SU (3) l ® SU (3)^, chiral transformation, 
the basis operators Qi transform as (27, 1)© (8, 1) 
for i = 1, . . . , 4, 9, 10, (8, 1) for i = 5, 6 and (8, 8) 
for i — 7, 8. To leading order in xPT the op- 
erators that transform according to these flavor 
representations are given byj^-jl] 

(8,1): A = (a^Eta^E)23, (2) 
B =(S+M + M+S)23 (3) 
(27,1): C =3(Sta,,S)23(Sta^S)ii 

+2(sta^s)i3(i]ta^i])2i (4) 

(8,8): D =st,Ei3, (5) 

where we consider (8,8) operator only at 0{p^) 
order due to the reason explained later. 

Using unknown parameters, a,b,c,d. . ., one 
may write the Qj's in terms of the xPT opera- 
tors according to 



Qi = aiA + b^B + CiC {i 
Q, = (i = 7,8) 



1, 



.,6,9,10) (6) 
(7) 



for the case of / = and 2 separately. 

For Qi in (^, 7r+ matrix elements are 

proportional to ai-\-hi~2ci, whereas tt^tt~ 
matrix elements to a.; — 2ci. To eliminate bi one 
has to introduce an operator 

Qsub = ("Is + rnd)sd - {nis - mdjs-fzd = 6subB (8) 

and construct Qi — aiQsub, where ai — bi/bsuh 
can be determined by —> matrix elements 
(See (ffl])). One then obtains the reduction for- 
mula, (pf) below, relating the K — > tttt and K tt 
amplitudes originally given by Bernard et. al^ . 

In (|^) there would appear too many unknown 
parameters to invalidate the reduction relation, 
if all (8,8) operators up to 0{p'^) order were 
included as in the case of other representation. 
To exclude contributions from these operators at 
O(p^), we take the chiral limit where the reduced 
relation is justified. In this limit, the matrix ele- 
ments for K — > TTTT are obtained in proportion to 
those for is:^ 7r|,|. 

It should be stressed that it is only in the chi- 
ral limit that one can obtain a complete set of 
reduction formulae at the lowest order of xPT. 



Table 1 

Simulation parameters. Scale a" 
string tension y/a = 440 MeV. 



^ is fixed from 



gauge action 


plaquette 


RG-improved 


P 


6.0 


2.6 


a-' 


2.0 GeV 


1.94 GcV 


lattice size 


16^ X 32 X 16 


DW height M 


1.8 


ruf 


0.01 - 0.06 in steps of 0.01 


^ji^config. 


111 


110 


"I5g 10 


2.96(34) McV 


0.283(42) MeV 



3. NUMERICAL SIMULATION 

Parameters of our numerical simulation are 
summarized in Table. |^. Calculations are car- 
ried out, assuming degenerate bare quark mass 
ruf — viu = rud = rus, on a 16"^ x 32 x 16 lat- 
tice at a coupling constant corresponding to a lat- 
tice spacing of « 2 GeV. With these parame- 
ters the anomalous quark mass m<^q representing 
residual chiral symmetry breaking is similar to u- 
d quark mass for the plaquette action and quite 
small for the RG-improved action. We set the 
scale from measurement of the string tension as- 
suming ^ = 440 MeVQ. 

To evaluate quark loops in some types of con- 
tractions, we employ the random noise method 
with the number of noises taken to be 2. This 
number is decided from a numerical test. 

We calculate the if — > tt matrix elements using 
the wall source and dividing by the normalization 
factor (7r+ IA4I 0) (0 IA4I K+). In this convention. 



xPT reduction formulae for 
take the form, 



1, 



.,6,9,10 



Q 



(I) 



(77+ 






(7r+ 


\A4\0) (0IA4I. 


K+) 



(9) 



and for i = 7,1 



Q 



(!) 



-V2U 



X TO 



Q 



M' 



,(10) 



IA4IO) {<d\Ai\K+)-' 

where vtlm is the meson mass on the lattice, while 
the experimental values are to be substituted in 
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(b) RG improved 
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Figure 1. Propagators relevant to Qg*^-* with 
nif — 0.03 for (a) plaquette and (b) RG- improved 
gauge actions. Fit ranges and jack knife errors are 
represented by horizontal lines. 



, m^r , JTiif . In (||) the subtraction term a^Qsub 
appears only for / = 0. The parameter ai is 
written as 









( 


3 \sj5d\ RO) 



(11) 



where the differentiation by rris is implemented 
by a double inversion of the Dirac operator. 

4. RESULTS 

For limitation of space we concentrate on the 
matrix elements of Q^^^ and Qg^"*, which represent 



the main contributions to e'/e. 

Examples of propagator ratios for Q^^'^ which 
appears on the right hand-side of (|^) are shown 
in Fig. 1^ for (a) plaquette and (b) RG-improved 
gauge actions. In each figure, the upper panel 
is for the K ^ tt matrix element and the lower 
for ae defined in (ITTI). Although both propaga- 
tors contain quark loops, which have been known 
to make signals worse, reasonable signals are ob- 
tained at the level of propagator ratios. This also 
justifies our choice of 2 random noises to evaluate 
the quark loops. 

A constant fitting, as shown in Fig. |l|, yields the 
matrix element, which we plot as a function of mj 
in Fig. ^. Circles and diamonds show the contri- 
bution of the if — » TT matrix element and the sub- 
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Figure 2. (5g°^ calculated by the xPT reduction 
formula (^. Original matrix element (circles), 
subtraction term (diamonds) and the total value 
(squares) are separately plotted. Lines are guide 
to eyes. 

traction term to the K — > tttt matrix element of 
Qg°' as given in the reduction formula . There 
is a severe cancellation between the two contribu- 
tions, leading to the total value of the K ^ tttt 
matrix element (squares) which is more than an 
order of magnitude smaller than the individual 
contributions. Nonetheless, the enlarged plots in 
Fig. ^(a) and (b) show that the total values have 
a reasonable signal. 

(2) 

There is no subtraction term for Qg ' . We ob- 
tain clear signals for this operator using the re- 
duction formula ( p^ ) as shown in Fig. |[ 

The remaining procedure is to make a linear 
fitting, as shown in Figs. ^ and ^, to take the 
chiral limit m/ — > where both of the reduc- 
tion formulae ^ and ( |To| ) are valid to estimate 
the physical matrix element, and to convert the 
values to continuum theory renormalized in the 
MS scheme with NDR. For the latter, we employ 
the renormalization factors calculated in pertur- 
bation theory at one- loop level with mean field 
improvement. 

Prior to this final step, we are currently increas- 
ing the statistics to reduce the errors further. We 
are also performing a new simulation with a larger 
lattice volume to investigate finite size effect in 
the matrix elements. 
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Figure S. K ^ tttt matrix element of Qg*^"* and linear extrapolation to the limit m/ — > for the plaquette 
and RG-improved gauge action. 



(a) plaquette 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 



(b) RG improved 



Figure 4. Counterpart of Fig. ^ for the operator Qg^' 

5. OUTLOOK 

We have presented our preliminary results for 
the K — !■ TTTT decay amplitudes based on the 
xPT reduction formulae and domain wall QCD. 
Matrix elements of reasonable statistical quality 
are obtained from about a hundred gauge config- 
urations in our quenched numerical simulation. 
Values for the matrix elements from plaquette 
and RG-improved gauge actions seems consistent 
within the errors. These results make us hopeful 
that the present approach yields precise informa- 
tion about the AI =1/2 rule and direct CP vio- 
lation in the standard model with more statistics 
and detailed analysis. 
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